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The event-by-event px— fluctuations due to the temperature variations are considered. The so- 
called measure is computed and shown to be a linear function of the temperature variance. The 
fluctuations of temperature naturally explain the data on $(p±) in proton-proton and central Pb-Pb 
collisions but independent measurements of the temperature fluctuations are needed to confirm the 
, explanation. Feasibility of such an event-by-event measurement is discussed. 
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Predictions of different models of heavy-ion collisions are often quite similar when averaged characteristics of the 
^ ' collisions are considered. Fluctuations are usually much more sensitive to the collision dynamics and consequently 
can be helpful in discriminating among the models. Since large acceptance detectors, which have recently become 
I common, make possible a detailed analysis of individual collisions the study of the event-by-event fluctuations appears 

■ to be a very promising field of high-energy heavy-ion physics, see for a review. 
The fluctuations in proton-proton and central Pb-Pb collisions at 158 GeV per nucleon have been recently 

^ measured Q on event-by-event basis. To eliminate trivial 'geometrical' fluctuations due to the impact parameter 

variation, the so-called measure Q has been used. $ is constructed is such a way that it is exactly the same 
for nucleon-nucleon (N-N) and nucleus-nucleus (A-A) collisions if the A-A collision is a simple superposition of the 
N-N interactions. In that case $ is independent of the centrality of an A-A collision. Moreover, $ equals zero 
when the inter-particlc correlations are entirely absent. The critical analysis of the $— measure can be found in [|,|. 
Q In the central Pb-Pb collisions the measured value of equals 4.6 ±1.5 MeV while the preliminary result for 

^ . proton-proton interactions in the same acceptance region is 5 ± 1 MeV ||]. Although the two values are close, the 
^ ■ mechanisms behind them seem to be different. It has been shown j^] that the correlations, which have short range in 
^ , the momentum space, like those due to the quantum statistics, are responsible for the positive value of ^(px) in the 
central Pb-Pb collisions. When these correlations are subtracted = 0.6 ± 1.0 MeV Q|. Our calculations have 

. , also demonstrated that the effect of Bose statistics of pions reduced by the hadron resonances fully explains the 

■ experimentally observed ^{p±) — 4.6 ± 1.5 in the central Pb-Pb collisions. In the p-p case, the situation seems to 
" ■ ' be opposite - the short range correlations provide a negligible contribution to ^{p±_) while the whole effect is due to 

the long range fluctuations. Thus, the data suggest that the dynamical long range correlations are reduced in the 
central Pb-Pb collisions (when compared to p-p) with the short range caused by the Bose statistics being amplified. 
The former feature is a natural consequence of the system's evolution towards the thermodynamic equilibrium. The 
amplification of the quantum statistics effect results from the increased particle population in the final state phase- 
space. Since various dynamical correlations contribute to ^(p±) the question arises what is the dynamical correlation 
in the nucleon-nucleon interactions which appears to be absent in the central nucleus-nucleus collisions. 

In the recent paper of the two of us the correlation, which couples the average p± to the event multiplicity 
N, has been studied. The correlation is convincingly evidenced in the p-p collisions The approximate analytical 
formula of ^{jpi_) as a function of the correlation strength has been derived and then the numerical simulation has 
been performed taking into account the finite detector's acceptance. The effect of the correlation (p^) vs. N has been 
shown to be very weak if the particles from a small acceptance region are studied. Consequently, the correlation is 
far too small to explain the preliminary experimental value of $(p_l) in the proton-proton collisions 

Our aim here is to discuss another possible mechanism responsible for the finite value of $(px) in p-p interactions 
[^. Namely, we analyze the effect of the temperature fiuctuations. Its role in shaping the particle spectra has been 
studied before pO| . Here, the temperature, or more generally the slope parameter of the p^— distribution, is assumed 
to vary from event to event. We compute ^(px) and find it to be a linear function of the temperature variance. As is 
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well known , the T— variance is directly related to the system's heat capacity. The idea to exploit the relationship 
to determine the specific heat of matter produced in nuclear collisions has been formulated in ||l2|,|l3|, see also [Of. 
Since the temperature fluctuations have been shown ||l5| to yield the so-called non-extensive Tsallis statistics 
(with a power-law instead of an exponential energy distribution) we express ^{pi_) by the non-extensivity parameter 
q. Further, we perform the numerical simulation of the p-p interactions with the effect of detector's finite acceptance 
taken into account. The temperature fluctuations are shown to explain in a natural way the data on ^{p±_) in proton- 
proton and central Pb-Pb collisions. Finally, we discuss how to perform an independent measurement of the genuine 
temperature fluctuations which have to be extracted from the statistical noise. 



II. ANALYTICAL CONSIDERATIONS 

Let us first introduce the $— measure. One defines a single-particle variable z = x — x with the overline denoting 
average over a single particle inclusive distribution. Here, we identify x with p± - the particle transverse momentum. 

The event variable Z, which is a multiparticle analog of z, is defined as Z '^iLii^i ~ where the summation 
runs over particles from a given event. By construction, (Z) = 0, where (...) represents averaging over events. Finally, 
the measure is defined in the following way 



Various fluctuations or correlations contribute to (|l|). Our aim here is to compute $(p±) when the temperature 
varies from event to event. P[T)iP±) denotes the single particle transverse momentum distribution in events with 
temperature T which is assumed to be independent of the event's multiplicity N. As discussed in [^^, the 
temperature can vary within the event but we discard such a possibility and assume that there is single temperature 
characterizing every event. We will return to this point in the concluding section. Then, the inclusive transverse 
momentum distribution reads 

/"CXD 

Pincl{p±) - / dT V{T) P^T){P±) , (2) 

Jo 

where V{T) describes the temperature fluctuations. Consequently, 

dT V{T) / dpx {p± - pifPiT) {P±) , (3) 
JO Jo 

with 
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PT= dTViT) dp^ p^ P^T){P±) ■ 
Jo Jo 



The A^— particle transverse momentum distribution in the events of multiplicity N is assumed to be the A^— product of 
P(t){P-l) weighed by the multiplicity and temperature distributions. Therefore, all inter-particle correlations, different 
than those due to the temperature variations, are entirely neglected. Then, one finds 
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Z^) =y,VN dTV{T) dpi P^T) (pi) ... / dpi P^T) {pf)(p\ + ...+pf- Npl) , 
Jo Jo Jo ^ ' 
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where Vj^ is the multiplicity distribution. Although the multi-particle distribution, Eq. (Q), may look as a simple 
product of the one-particle distributions, the particle distributions are not independent from each other due to the 
integration over T . 

In our further calculation we choose P(^T)iPi-) in ttic form suggested by the thermal model i.e. 



P(T)iP±) --P±exp 



y/m'^ + p\ 
T 



(5) 



where m is the particle mass. If the transverse collective flow is taken into account T should be understood as an 
effective temperature or simply a slope parameter controlled by the actual freeze-out temperature and the collective 
flow velocity. 
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In the limit m = the pj^— distribution (^) acquires a simple exponential form and one easily computes and 
(Z^) given by Eq. (||) and (^), respectively. Then, one gets 

i2 = 2(T2) +4((r2) - (T>2) 

which gives 



(iV> 



(T) 



(6) 



when the T— fluctuations are sufflciently small i.e. {N){T)^ 3> {N^){{T'^) — {T)"^). For the poissonian multiplicity 
distribution the formula simplifies to 



(7) 



III. INTERPRETATION OF -l- 

If T in Eq. corresponds to genuine temperature, not a slope parameter, the formula gets a nice interpretation 
due to the well-known thermodynamical relation |Tl| ] which has been discussed in the context of nuclear collisions in 
Namely, 

l_ ^ (T^) - (r)2 
a (T)2 ' 

where is the system's heat capacity. We note that C^, as an extensive thermodynamical parameter, is proportional 
to (iVtot) which is the average number of all particles - charged and neutral - in the system at freeze-out, not to the 
average number of the observed particles {N) which enters Eq. (^. Therefore, the formula can be rewritten using 
the relation (||) as 

with Cv = Cv/{Ntot) being the specific heat of hadronic matter at freeze-out. For a system of massless non-interacting 
bosons with vanishing chemical potential c„ = 27r"'/15C(3) = 10.8. We note that c„ is independent of the particle's 
internal degrees of freedom. 

A comment is in order here. It has been shown by one of us ^ that ^{p±) vanishes in the ideal classical gas, where 
the particles are exactly independent from each other. On the other hand, the formula (||) tells us that > in 

such a gas. However, there is no conflict between the two results. The computation from was performed at fixed 
temperature while the thermodynamical identity (|8|) states that there are temperature fiuctuations in any system with 
finite heat capacity. In fact, these fluctuations are usually very small, because the temperature variance is inversely 
proportional to the number of particles. Consequently, the variations of temperature are neglected in most cases. 
However, ^{p±) appears to be very sensitive to the temperature fluctuations and the two results differ qualitatively. 

In the recent paper of one of us , the so-called non-extensive Tsallis statistics has been shown to naturally 
emerge when a system experiences temperature fiuctuations. Specifically, it has been argued [ p^ that l/T often varies 
according to the gamma distribution. Then, 

^(T) = ^(-) exp(--), (10) 

with the parameters A and a related to the moments of 1 /T as 
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Substituting Eqs. (||,|To|) into one gets an inclusive distribution in the Tsallis statistics form ||Tf 

PhicliPA 



L -to 



(11) 



where Tq = (l/T) ^ and q — (A + 1)/A is the non-extensivity parameter |l6[ related to the temperature fluctuations 
as lUl 



(12) 



The second approximate equality holds for sufficiently small fluctuations. As known | [T^ , the distribution (^Tj) tends 
to with T = To when q ^ 1. Using the relation ([1^), the formula can be rewritten in yet another form, 



(13) 



which relates $ to the non-extensivity parameter q. 

The relationship between the measure and Tsallis parameter q has been earlier considered in a different context 
1^. Namely, the authors have studied how the g— statistics modifies the usual Bose-Einstein correlations discussed 
^. Then, $ has been found to decrease, not increase as in Eq. with q. 
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FIG. 1. $(pt) as a function of the temperature variance for three values of the acceptance probability p. The temperature 
varies according to the gaussian distribution while the multiplicity is controlled by the lognormal one. 



IV. NUMERICAL SIMULATION 



In this section we present the results of our Monte Carlo simulation of p-p collisions. The single-particle 
px— distribution is still given by Eq. (|^). The mass equals now that of a charged pion because all particles in 
our simulation are treated as charged pious. We have considered two temperature distributions: the gamma-like form 
( p^ ) and a gaussian distribution (cut off at T < 0). The lognormal distribution of multiplicity of negative particles 
has been shown to fit the p-p data very well in the broad range of the collision energies 18 1. We have used the 
parametrization given in |18|| and assumed that the numbers of positive and negative pious are equal to each other 
in every event. The assumption is certainly reasonable in the central rapidity domain. To check whether the results 
are sensitive to the form of the multiplicity distribution we have also performed a simulation with the poissonian 
distribution of negative particles. As in the case of the lognormal distribution, the multiplicity of charged particles 
has been simply assumed to be two times bigger than that of negative ones. The average charged particle multiplicity 
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and temperature have been taken as in our previous paper ^ i.e. (T) = 167 MeV and (N) = 6.56. These values 
correspond to the proton-proton collisions at 205 GeV. 
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FIG. 2. ^{pt) as a function of the temperature variance for three values of the acceptance probability p. The inverse 
temperature varies according to the gamma distribution while the multiplicity is controlled by the Poisson distribution. 

Due to the particle registration inefficiency and finite detector's coverage of the final state phase-space, only a 
fraction of the produced particles is usually observed in the experimental studies. Our Monte Carlo simulation takes 
into account the two effects in such a way that each generated particle - positive or negative pion - is registered 
with probability p and rejected with (1 — p). The detector's acceptance usually covers a given rapidity window but 
within our model the T— fluctuations and p^— distribution are rapidity independent. Therefore, there is no difference 
between a particle being lost due to the limited acceptance and one lost due to the tracking inefficiency. 




<N> 

FIG. 3. ${pt) as a function of the average number of observed particles for three values of the temperature variance. The 
temperature varies according to the gaussian distribution while the multiplicity is controlled by the lognormal one. 

The results of our p-p simulation arc shown in Figs. 1-4. Those in Figs. 1 and 3 have been found with the gaussian 
distribution of temperature and the lognormal multiplicity distribution. The results from Figs. 2 and 4 correspond 
to the gamma and Poisson distributions, respectively. When the gamma distribution is used the T— variance divided 
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by (T)^ is denoted by g — 1, in agreement with Eq. In the case of gaussian distribution the same quantity is 

written as cP' j (T)^. One sees in Figs. 1 and 2, that $(pt) grows hnearly with the temperature variance, exactly as in 
Eq. (^. As can also be seen, the gaussian and gamma distributions yield very similar results. The two observations 
mean that the effect of finite pion mass is small and that <I>, as in the m = case, is simply a linear function of the 
second moment of T . 
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FIG. 4. "1>(pt) as a function of the average number of observed particles for three values of the temperature variance. The 
inverse temperature varies according to the gamma distribution while the multiplicity is controlled by the Poisson distribution. 

Instead of the acceptance parameter p one can use the average multiplicity of the observed particles (N) to charac- 
terize the acceptance. In Figs. 3 and 4 we present $(j5t) as a function of (N). The growth of $(j3t) with (N) is seen 
to be almost linear. 



V. COMPARISON WITH THE EXPERIMENTAL DATA 

The preliminary experimental value of ^{pr) in proton-proton collisions is, as already mentioned, 5 ± 1 MeV [1). 
The measurement has been performed in the transverse momentum and pion rapidity intervals (0.005, 1.5) GeV and 
(4.0, 5.5), respectively. Only about 20% of all produced particles have been observed. According to our simulation 
one needs (t'^/{T)'^ = g - 1 ^ 0.015 ± 0.003, which corresponds to y^{T^) - (r)2 = 20 ± 3 MeV at (T) = 167 MeV, to 
reproduce the experimental result. 

Let us now calculate the specific heat of the hadronic matter produced in the proton-proton interactions from the 
obtained value of a/ (T^) — (T)-^. For that we identify (iVtot) from Eq. (^ with the total number of pions. The pions 
include those 'hidden' in hadron resonances. We count each p as two pions, each uj as three, etc. Then, Cy from (|^) 
is the heat capacity per pion. We estimate (iVtot) as 5 • 1.5 • (N), where the factor 5 is due to the acceptance and 1.5 
to include the neutral particles. Then, one finds from Eq. (0) the heat capacity per pion Cy — 6 zt 2. This number 
is significantly smaller than the previously mentioned specific heat of massless non-interacting bosons which equals 
10.8. In fact, the discrepancy is even worse because a more realistic model of strongly interacting matter which takes 
into account numerous resonances and, obviously, finite hadron masses gives the heat capacity per pion exceeding 20 
[ p^ . However, the very applicability of the thermodynamical model to the p-p collisions is far from obvious. 

As discussed in the Introduction, ^{p±), which is measured in the central Pb-Pb collisions, equals 4.6 ± 1.5 MeV 
This value includes the short range (Bose-Einstein) correlations. When those correlations are excluded $(p_l) = 
0.6 ±1.0 MeV ||]. Since the effect of quantum statistics is not taken into account in our simulation, this is the latter 
experimental value of $(pj_) which should be compared with our calculation. The observed average multiplicity has 
been 270, i.e., as in the p-p interactions, about 20% of all produced charged particles ||]. We first identify the system 
freeze-out temperature with the slope parameter deduced from the pion transverse momentum distribution (T) = 180 
MeV |l9). Then, $(pi) = 0.6 ± 1.0 MeV yields via Eq. (0) y/{T'^) - (T)2 = 0.5 ± 0.4 MeV. Let us stress here that 



our numerical simulation fully confirms the reliability of the analytical formula (|7|). The temperature is significantly 
reduced if the transverse hydrodynamic expansion is taken into account. The freeze-out temperature obtained by 
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means of the simultaneous analysis of the single particle spectra and the Bose-Einstein (HBT) correlations is about 
120 MeV This temperature combined with <^{p±) = 0.6 ± 1.0 MeV and (N) = 270 gives ^y{T^) - (T)2 = 

0.4 ± 0.4 McV. It has been argued in that the transverse collective flow, which significantly modifies the observed 
temperature, does not contribute much to the slope parameter fluctuations. Therefore, the temperature variance is 
expected to be dominated by the genuine temperature fluctuations. 

Calculating, as in the p-p case, the heat capacity per pion from the temperature dispersion, one finds that (T) = 180 
MeV and y/{T'^) - (T)2 = 0.5 ± 0.4 MeV give Cy = 60 ± 100 while (T) = 120 MeV and y/{T^) - (T)2 = 0.4 ± 0.4 
MeV correspond to Ct, = 40 ± 70. Unfortunately, the errors are so large that no conclusion can be drawn. 



VI. MEASUREMENT OF T FLUCTUATIONS 



In the previous section we have shown that the temperature fiuctuations naturally explain the p-p and Pb-Pb data. 
Here, let us briefiy consider how to observe independently the event-by-event temperature fluctuations. We discuss 
a straightforward method proposed in |po[| . However, other procedures, in particular the so-called sub-event method 
developed in ||2l]] , might be more efficient. The temperature variance can be found measuring the event's average 
transverse mass defined as 

1 ^ 



where N denotes the event's multiplicity and m^^ ~ ^ + p^^ is the transverse mass of i— th particle. If the single 
particle p_L— distribution is of the form (||) the m^— distribution reads 



-P(T)(™±) ~ ™±exp[- — ] (14) 



and /i^ is related to T in the following way 

2^2 2Tm + 

r + m 



Then, the event's temperature is expressed through as 
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T = - [y4TO/Xi - 4m2 + + /i_L - 2m] . (15) 

Thus, measuring /^j^ on the event-by-event basis one can get the temperature variance V{T) = (T^) — (T)^. However, 
the statistical fluctuations due to the finite event multiplicity have to be subtracted. The point is that when the 
genuine temperature does not fluctuate at all, the observed temperature does vary because the number of registered 
particles is not infinite. 

When the genuine temperature is fixed and the particles are independent form each other the variance of /i^ is 
fully determined by the statistical fluctuations. In the events of multiplicity N it equals 

where Vi(tox) computed with the m^— distribution ( p^ is 

, + QT'^m + ■^Tm'^ + /2T'^ + 2Tm + \^ 

Vi{m±) = — — . (16) 

1 + m \ 1 + m I 

The T— variance is found from the /i^— variance in the usual way [ p2| i.e. 

^^(r)-f:^)V.(Mx), (17) 



where the derivative is computed at T = (T) and /i^ = Since T is not a linear function of Eq. (^^ holds for 

sufficiently small 14 (m±) |2^- This in turn demands that ^ 1. Using Eq. (p^), one finally finds the contribution 
of statistical fluctuations to the temperature variance as 
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2m + 



N 



(18) 



where Vi {m± ) is given by Eq. (|l|) with T replaced by (T). We note that Vs{T) = {T)'^/N when m = 0. The variance 
Vs{T) should be subtracted from the observed T— variance V{T) to eliminate the statistical fluctuations. 
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FIG. 5. The observed T— variance as a function of multiplicity. 
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FIG. 6. The T— variance with the subtracted statistical fluctuations as a function of the multiplicity. 



We have performed a simple simulation to see how well the subtraction procedure works. For this purpose we have 
generated the events with fixed multiplicity N and temperature fluctuating from event to event. The single particle 
p_L— distribution has still been of the form (|^) and the temperature varied according to the gaussian distribution with 
(T) = 167 MeV. The event temperature has been determined in the 'experimental' way described above. Finally, 
the temperature variance has been found and the statistical contribution has been subtracted. In Figs. 5 and 6 we 
show V{T) and V{T) — Vs{T) as a function of the event multiplicity. One can see that in spite of large values of 
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V(T) and Vs(T) the input temperature variance = 100 MeV^ is well reproduced for N as small as 10. The result, 
although expected, is not entirely trivial. The expression ( |l8|) assumes that the particles are independent from each 
other while in our simulation the particles are correlated because of the temperature fluctuations. We conclude that 
the measurement of T— fluctuations seems to be a feasible task even in the case of relatively low multiplicity collisions. 



VII. FINAL REMARKS 

We have shown that $(j3±) observed in proton-proton collisions Q can be understood as an effect of temperature 
fluctuations with (T^) — (T)^ = 20 MeV. While the result needs to be confirmed by independent T— variance 
measurements let us mention here an interesting observation from p3t . It has been found there that the transverse 
mass spectrum of ir^ from p-p collisions at ^/s — 30 GeV decreases as pj^^ over 10 orders of magnitude with 
P = 9.6. Within the thermal model, such a behavior naturally appears due to the temperature fluctuations p5[. 



Then, the exponent P — 9.6 gives the non-extensivity parameter in Eq. (11) equal q — 1.10 which translates into 
-\/ (T^) — (T)2 — 53 MeV at (T) = 167 MeV. The two values of the temperature dispersion extracted from the p-p 
data, 20 MeV and 53 MeV, have been found in different ways. One easily shows that $ is mostly sensitive to the 
event-by-event temperature fluctuations while the inclusive distribution is shaped both by the temperature variation 
from event to event and within the event. Therefore, the T— dispersion found from the inclusive spectrum is expected 
to be larger than that from 

^{p±) observed in the central Pb-Pb collisions is very small when the short range correlations are subtracted [1). 
This small value corresponds to the temperature dispersion below 1 MeV. The difference of the dispersions found in 
p-p and central Pb-Pb collisions is not surprising. The hadronic system from nuclear interactions is not only larger - 
$ as an intensive quantity is not directly dependent on the system size - but at freeze-out it is expected to be much 
closer to the thermodynamic equilibrium. Consequently, the temperature fluctuations should be smaller. 

We conclude our considerations as follows. The evcnt-by-event fluctuations of temperature are a possible mechanism 
determining the value of ^{p±) in p-p collisions. The smallness of the contribution of the long range correlations to 
^{pi_) in the central Pb-Pb collisions is then also naturally explained. One needs an independent measurement of 
the temperature variance to confirm the explanation. Since the sizable fluctuations due to the finite statistics are, in 
principle, under control, such a measurement seems to be feasible even in the relatively low multiplicity interactions. 
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